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The effect of convection on pulsational stability 
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Abstract 

A review on the current state of mode physics in classical pulsators is presented. 
Two, currently in use, time-dependent convection models are compared and 
their applications on mode stability are discussed with particular emphasis on 
the location of the Delta Scuti instability strip. 

Introduction 

Stars with relatively low surface temperatures show distinctive envelope convec- 
tion zones which affect mode stability. Among the first problems of this nature 
was the modelling of the red edge of the classical instability strip (IS) in the 
Hertzsprung-Russell (H-R) diagram. The first pulsation calculations of classi- 
cal pulsators without any pulsation-convection modelling predicted red edges 
which were much too cool and which were at best only neutrally stable. What 
follows were several attempts to bring the theoretically predicted location of 
the red edge in better agreement with the observed location by using time- 
dependent convection models in the pulsation analyses (Dupree 1977; Baker &i 
Gough 1979; Gonzi 1982; Stellingwerf 1984). More recently several authors, 
e.g. Bono et al. (1995, 1999), Houdek (1997, 2000), Xiong & Deng (2001, 
2007), Dupret et al. (2005) were successful to model the red edge of the clas- 
sical IS. These authors report, however, that different physical mechanisms are 
responsible for the return to stability. For example. Bono et al. (1995) and 
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Figure 1: Sketch of an overturning hexagonal 
(dashed lines) convective cell. Near the centre 
the gas raises from the hot bottom to the cooler 
top (surface) where it moves nearly horizontally 
towards the edges, thereby loosing heat. The 
cooled gas then descends along the edges to close 
the circular flow. Arrows indicate the direction of 
the flow pattern. 
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Table 1: Summary of time-dependent convection model differences. 



Balance between buoyancy & 
turbulent drag (Unno 1967, 1977) 

- acceleration terms of convective 
fluctuations w, T' neglected 

- nonlinear terms approximated 
by spatial gradients (xl/t 

- dip' neglected in momentum equ. 

- characteristic eddy lifetime: 

T ~ i/2w 

- variation (.i = 51/1 (Unno 1967): 
WT < 1 : £i ^ Hi 

LOT > 1 : ^1 ~ ri 
or (Unno 1977): 

£i - (1 + ia;V2)-i(iJi - iwT/9i/3) 
{H is pressure scale height) 

- turbulent pressure pt neglected 
in hydrostatic support equation 



Kinetic theory of accelerating 
eddies (Gough 1965, 1977a) 

- acceleration terms included: w, 
T' evolve with growth rate a 

- nonlinear terms are neglected 
during eddy growth 

- dip' included in Eq. (1) 

- r = 2/(7 determined stochas- 
tically from parametrized shear 
instability 

- variation of mixing length 
according to rapid distortion 
theory (Townsend 1976), i.e. 
variation also of eddy shape 



- Pt = pww included in mean 
equ. for hydrostatic support 



Dupret et al. (2005) report that it is mainly the convective heat flux, Xiong &l 
Deng (2001) the turbulent viscosity, and Baker &i Gough (1979) and Houdek 
(2000) predominantly the momentum flux (turbulent pressure pt) that stabilizes 
the pulsation modes at the red edge. 

Time-dependent convection models 

The authors mentioned in the previous section used different implementations 
for modelling the interaction of the turbulent velocity field with the pulsation. In 
the past various time-dependent convection models w/ere proposed, for example, 
by Schatzman (1956), Gough (1965, 1977a), Unno (1967, 1977), Xiong (1977, 
1989), Stellingwerf (1982), KuhfuB (1986), Canute (1992), Gabriel (1996), 
Grigahcene et al. (2005). Here I shall briefly review and compare the basic 
concepts of two, currently in use, convection models. The first model is that 
by Gough (1977a, b), which has been used, for example, by Baker &i Gough 
(1979), Balmforth (1992) and by Houdek (2000). The second model is that by 
Unno (1967, 1977), upon which the generalized models by Gabriel (1996) and 
Grigahcene et al. (2005) are based, with applications by Dupret et al. (2005). 
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Nearly all of the time-dependent convection models assume the Boussinesq 
approximation to the equations of motion. The Boussinesq approximation relies 
on the fact that the height of the fluid layer is small compared with the density 
scale height. It is based on a careful scaling argument and an expansion in small 
parameters (Spiegel Veronis 1960; Gough 1969). The fluctuating convection 
equations for an inviscid Boussinesq fluid in a static plane-parallel atmosphere 
are 

dtUi + {ujdjUi - UjdjUi) = -p~^dip' + gaT' Si3 , (1) 
dtT' + [ujdjT' -^^^) = 0w-{pc^)-^diF!, (2) 

supplemented by the continuity equation for an incompressible gas, djUj = 0, 
where u = {u,v,iu) is the turbulent velocity field, p is density, p is gas pres- 
sure, g is the acceleration due to gravity, T is temperature, Cp is the specific 
heat at constant pressure, a = —{dlnp/dlnT)p / T, Fi is the radiative heat 
flux, P is the superadiabatic temperature gradient and 5ij is the Kronecker 
delta. Primes (') indicate Eulerian fluctuations and overbars horizontal aver- 
ages. These are the starting equations for the two physical pictures describing 
the motion of an overturning convective eddy, illustrated in Fig. 1. 
In the first physical picture, adopted by Unno (1967), the turbulent element, 
with a characteristic vertical length £, evolves out of some chaotic state and 
achieves steady motion very quickly. The fluid element maintains exact balance 
between buoyancy force and turbulent drag by continuous exchange of momen- 
tum with other elements and its surroundings. Thus the acceleration terms dfUi 
and dtT' are neglected and the nonlinear advection terms provide dissipation (of 
kinetic energy) that balances the driving terms. The nonlinear advection terms 
are approximated by UjdjUi — UjdjUi ~ 2w'^ /£ and UjdjT' — UjdjT' IwT' jl. 
This leads to two nonlinear equations which need to be solved numerically to- 
gether with the mean equations of the stellar structure. 

The second physical picture, which was generalized by Gough (1965, 1977a, b) 
to the time-dependent case, interprets the turbulent flow by indirect analogy 
with kinetic gas theory. The motion is not steady and one imagines the con- 
vective element to accelerate from rest followed by an instantaneous breakup 
after the element's lifetime. Thus the nonlinear advection terms are neglected 
in the convective fluctuation equations (l)-(2) but are taken to be responsi- 
ble for the creation and destruction of the convective eddies (Gough 1977a, b). 
By retaining only the acceleration terms the equations become linear with an- 
alytical solutions w oc exp((jf) and T' oc exp{at} subject to proper periodic 
spatial boundary conditions, where t is time and ^R{(7) is the linear convec- 
tive growth rate. The mixing length £ enters in the calculation of the eddy's 
survival probability, which is proportional to the eddy's internal shear (rms vor- 
ticity), for determining the convective heat and momentum fluxes. Although 
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Figure 2: Mode stability of an 1.7 Delta Scuti star computed with Cough's 
(1977a, b) convection model. Left: Stability coefficient ri — cui/uj^ as a function 
of surface temperature Tctt across the IS. Results are shown for the fundamental ra- 
dial mode {n — 1) and for two values of the mixing-length parameter a. Positive rj 
values indicate mode instability. Right: Integrated work integral W as a function of 
the depth co-ordinate log(r) for a model lying just outside the cool edge of the IS 
{Taff = 6813 K). Results are plotted in units of r; and for a = 2.0. Contributions to W 
(solid curve) arising from the gas pressure perturbation, Wg (dashed curve), and the 
turbulent pressure fluctuations, Wt (dot-dashed curve), are indicated {W — W^+Wt). 
The dotted curve is the ratio of the convective to the total heat flux Fc/F. Ionization 
zones of H and He (5% to 95% ionization) are indicated (from Houdek 2000). 

the two physical pictures give the same result in a static envelope, the results 
for the fluctuating turbulent fluxes in a pulsating star are very different (Gough 
1977a). The main differences between Unno's and Cough's convection model 
are summarized in Table 1. 

Application on mode stability in 6 Scuti stars 

Fig. 2 displays the mode stability of an evolving 1.7 M0 Delta Scuti star crossing 
the IS. The results were computed with the time-dependent, nonlocal convec- 
tion model by Gough (1977a, b). As demonstrated in the right panel of Fig. 2, 
the dominating damping term to the work integral W for a star located near 
the red edge is the contribution from the turbulent pressure fluctuations Wt- 

Gabriel (1996) and more recently Grigahcene et al. (2005) generalized 
Unno's time-dependent convection model for stability computations of nonradial 
oscillation modes. They included in their mean thermal energy equation the 
viscous dissipation of turbulent kinetic energy, e, as an additional heat source. 
The dissipation of turbulent kinetic energy is introduced in the conservation 
equation for the turbulent kinetic energy K := UiUi/2 (e.g. Tennekes & Lumley 
1972, §3.4; Canute 1992; Houdek L Cough 1999): 



DtK + dj{Kuj + p^'^p'uj) - i^dfK = -u-uJdjUt + gau^T' - e , (3) 




Figure 3: Stability computations of Delta Scuti stars which include the viscous dissi- 
pation rate e of turbulent kinetic energy according to Grigahcene et al. (2005). Left: 
Blue and red edges of the IS superposed on evolutionary tracks on the theorists H-R 
diagram. The locations of the edges, labelled p„B and PnR, are indicated for radial 
modes with orders 1 < n < 7. Results by Houdek (2000, a — 2.0, see Fig. 2) and 
Xiong & Deng (2001) for the gravest p modes are plotted as the filled and open circles 
respectively. Right: Integrated work integral W as a function of the depth co-ordinate 
log(T) for a stable n = 3 radial mode of a 1.8 Mq star (see 'star' symbol in the left 
panel). Contributions to W arising from the radiative flux, W^r, the convective flux, 
Wc, the turbulent pressure fluctuations, Wt [Wr, + Wc + Wt), and from the pertur- 
bation of the turbulent kinetic energy dissipation, W^e {Wr + Wc + Wt), are indicated 
(adapted from Dupret et al. 2005). 

where is the material derivative, Ui is the average (oscillation) velocity, 
i.e. the total velocity Ui = Ui + ui, and v is the constant kinematic viscosity 
(in the limit of high Reynolds numbers the molecular transport term can be 
neglected). The first and second term on the right of Eq. (3) are the shear 
and buoyant productions of turbulent kinetic energy, whereas the last term 
e = ^{djUi + diUj)"^ /2 is the viscous dissipation of turbulent kinetic energy 
into heat. This term is also present in the mean thermal energy equation, but 
with opposite sign. The linearized perturbed mean thermal energy equation for 
a star pulsating radially with complex angular frequency uj ~ ojr + ioji can then 
be written, in the absence of nuclear reactions, as {'S' denotes a Lagrangian 
fluctuation and I omit overbars in the mean quantities): 

dSL/dm = -iujCpT{ST/T - V^dSp/p) + 5e , (4) 

where m is the radial mass co-ordinate, Vad = {d\nT/d\np)s and L is the to- 
tal (radiative and convective) luminosity. Grigahcene et al. (2005) evaluated e 
from a turbulent kinetic energy equation which was derived without the assump- 
tion of the Boussinesq approximation. Furthermore it is not obvious whether 
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the dominant buoyancy production term, gaUjT' (see Eq. 3), was included in 
their turbulent kinetic energy equation and so in their expression for e. 
Dupret et al. (2005) applied the convection model of Grigahcene et al. (2005) 
to Delta Scuti and 7 Doradus stars and reported well defined red edges. The 
results of their stability analysis for Delta Scuti stars are depicted in Fig. 3. 
The left panel compares the location of the red edge with results reported by 
Houdek (2000, see also Fig. 2) and Xiong & Deng (2001). The right panel of 
Fig. 3 displays the individual contributions to the accumulated work integral W 
for a star located near the red edge of the n = 3 mode (indicated by the 'star' 
symbol in the left panel). It demonstrates the near cancellation effect between 
the contributions of the turbulent kinetic energy dissipation , W^, and turbulent 
pressure, W^, making the contribution from the fluctuating convective heat flux, 
Wc, the dominating damping term. The near cancellation effect between TI7 
and Wt was demonstrated first by Ledoux &l Walraven (1958, §65) (see also 
Gabriel 1996) by writing the sum of both work integrals as: 

r-M 

We + Wt = 37r/2 / (5/3 - 73)3(<5ptMP"' dm , (5) 

Jmh 

where M is the stellar mass, mb is the enclosed mass at the bottom of the 
envelope and 73 = l-|-(91nT/31np)s (s is specific entropy) is the third adiabatic 
exponent. Except in ionization zones 73 ~ 5/3 and consequently W^ + Wt — 0. 

The convection model by Xiong (1977, 1989) uses transport equations for 
the second-order moments of the convective fluctuations. In the transport 
equation for the turbulent kinetic energy Xiong adopts the approximation by 
Hinze (1975) for the turbulent dissipation rate, i.e. e = 2xk{uiUip'^ /3p'^)^^'^ , 
where x = 0-45 is the Heisenberg eddy coupling coefficient and k oc is 
the wavenumber of the energy-containing eddies. However, Xiong does not 
provide a work integral for e (neither does Unno et al. 1989, §26,30) but 
includes the viscous damping effect of the small-scale turbulence in his model. 
The convection models considered here describe only the largest, most energy- 
containing eddies and ignore the dynamics of the small-scale eddies lying further 
down the turbulent cascade. Small-scale turbulence does, however, contribute 
directly to the turbulent fluxes and, under the assumption that they evolve 
isotropically, they generate an effective viscosity vt which is felt by a particular 
pulsation mode as an additional damping effect. The turbulent viscosity can 
be estimated as (e.g. Cough 1977b; Unno et al. 1989, §20) ia ^ XiwwY^'^l, 
where A is a parameter of order unity. The associated work integral Wi, can be 
written in Cartesian co-ordinates as (Ledoux & Walraven 1958, §63) 



Wi, = -27r cJr / H 

I nib 



"I 



Mr 1 

eii-3(V-0' 



dm , (6) 
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Figure 4: Accumulated work integral W as a function of the depth co-ordinate log(p). 
Results are shown for the n = 1 radial mode of a Delta Scuti star located inside 
the IS (left panel) and outside the red edge of the IS (right panel). The stability 
calculations include viscous dissipation by the small-scale turbulence (Xiong 1989; see 
Eq. 6). Contributions to W (solid curve) arising from the fluctuating gas pressure, 
Wg (dashed curve), the turbulent pressure perturbations, Wt (long-dashed curve), and 
from the turbulent viscosity, Wv (dotted curve), are indicated [W — Wg + Wt + Wv)- 
The ionization zones of H and He are marked (adapted from Xiong &i Deng 2007). 



where Cij = {dj£,i + di£,j)/2 and ^ is the displacement eigenfunction. Xiong &i 
Deng (2001, 2007) modelled successfully the IS of Delta Scuti and red giant 
stars and found the dominating damping effect to be the turbulent viscosity (Eq. 
6). This is illustrated in Fig. 4 for two Delta Scuti stars: one is located inside 
the IS (left panel), the other outside the cool edge of the IS (right panel). The 
contribution from the small-scale turbulence was also the dominant damping 
effect in the stability calculations by Xiong et al. (2000) of radial p modes in 
the Sun, although the authors still found unstable modes with orders between 
11 < n < 23. The importance of the turbulent damping was reported first 
by Goldreich & Keeley (1977) and later by Goldreich & Kumar (1991), who 
found all solar modes to be stable only if turbulent damping was included in 
their stability computations. In contrast, Balmforth (1992), who adopted the 
convection model of Gough (1977a, b), found all solar p modes to be stable 
due mainly to the damping of the turbulent pressure perturbations, Wt, and 
reported that viscous damping, W,y, is about one order of magnitude smaller 
than the contribution of Wt- Turbulent viscosity (Eq. 6) leads always to mode 
damping, where as the perturbation of the turbulent kinetic energy dissipation, 
Se (see Eq. 4), can contribute to both damping and driving of the pulsations 
(Gabriel 1996). The driving effect of 6e was shown by Dupret et al. (2005) for 
a 7 Doradus star. 
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Summary 

We discussed three different mode stability calculations of Delta Scuti stars 
which successfully reproduced the red edge of the IS. Each of these computa- 
tions adopted a different time-dependent convection description. The results 
were discussed by comparing work integrals. All convection descriptions in- 
clude, although in different ways, the perturbations of the turbulent fluxes. 
Gough (1977a), Xiong (1977, 1989), and Unnoetal. (1989) did not include the 
contribution Wc to the work integral because in the Boussinesq approximation 
(Spiegel &i Veronis 1960) the viscous dissipation is neglected in the thermal 
energy equation. In practise, however, this term may be important. Grigahcene 
et al. (2005) included Wc but ignored the damping contribution of the small- 
scale turbulence W„, which was found by Xiong Deng (2001, 2007) to be the 
dominating damping term. The small-scale damping effect was also ignored in 
the calculations by Houdek (2000). A more detailed comparison of the convec- 
tion descriptions has not yet been made but Houdek & Dupret have begun to 
address this problem. 
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Discussion 

Christensen-Dalsgaard: How does the mixing length affect the red edge of the 7 Dor 
instability strip? 

Houdek: The location of the red edge is predominantly determined by radiative 

damping which gradually dominates over the driving effect of the so-called convective 
flux blocking mechanism (Dupret et al. 2005). A change in the mixing length will 
not only affect the depth of the envelope convection zone but also the characteristic 
time scale of the convection and consequently the stability of g modes with different 
pulsation periods. A calibration of the mixing length to match the observed location 
of the 7 Dor instability strip will also calibrate the depth of the convection zone at a 
given surface temperature. 



